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Abstract 

A numerical modeling of the radiation emitted from a Luminescent dye embedded in a finite 
one-dimensional photonic crystal is presented. The Photonic Band Structure and the Photonic 
Density of States are derived using classical electromagnetic approach, and the Finite Difference 
Time-Domain formalism is used to calculate the electromagnetic field distribution. It is found 
that the periodic modulation provides an effective way to control the Spontaneous Emission under 
certain circumstances. We find the conditions where a large amount of light can be enhanced on the 
vicinity of a photonic band edge due to the presence of a high density of states. This phenomena 
opens the possibility to design new Lasers sources. 

PACS numbers: Valid PACS appear here 
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It is a long known fact that the Spontaneous Emission (SE) is not an immutable property 
of the source but can be altered by the environment in which the atom is situated [? ] . The 
SE in a periodic medium such as a Photonic Crystal (PC) leads to many effects not present 
on unbounded media [? ]. In particular, the modification of the radiation of a Luminescent 
dye embedded in a PC has been verified in many experiments. [? ? ? ] This phenomenon 
has considerable consequences in both, Science and Technology [? ] . In Science, because SE 
is a Quantum phenomenon effect that can be now study not only in relation to the existence 
of Photonic Band Gaps (PBG) -frequency ranges in which no propagating modes exist- but 
also within the interaction of more sophisticated phenomenas such as nonlinearities [? ? ], 
cavity interactions [? ], etc. For Technology, because a monochromatic enhancement of the 
emitted light opens the possibility to design new Laser sources [? ]. 

The SE have been explained by Classical and Quantum Electrodynamic approaches. 
[? ]. The Classical Electrodynamics approach study the SE in terms of the deformation 
of the radiation field of a dipole due to the multiple interferences of a periodic infinite 
media. [? ]. The Quantum Electrodynamics describes the SE as a formahsm based on 
the restriction of discrete electromagnetic field values when modeling it, and at first-order 
perturbation theory the emission rate can be calculated using Fermi's golden rule [? ]. In 
both approximations the formulation is based on the unrealistic consideration of an infinite 
crystal. A practical reahzation of SE on PC is inevitably of finite size, however. On the other 
hand, the consideration of an dipole as the source of radiation is also unrealistic. Dipolar 
radiation produces radiation with an intrinsic spatial distribution that does not occur on 
Luminescence experiments. [? ] The luminescent dye as a source of electromagnetic field is 
more likely to a Gaussian source with a central field frequency than the radiation field of a 
dipole. [? ] One has to note that that all kind of approximations. Classical or from Quantum 
theory have to establish its connection with the experimental conditions. In this case, the 
more suitable experiment to understand the Luminescence in presence of PBG materials is 
to study the SE modification of a Gaussian source immersed in a finite PC. 

The PBG materials are periodic structures which do not allow propagation of photons 
over a finite band of frequencies [? ]. A onc-dimcnsional (ID) PBG structure is the simplest 
PC and usually is called Bragg Reflector when restricted to a quarter- wave stack. [? ] ID- 
PBG structures have interesting effects as anomalous group velocity [? ] and dispersion [? 
] , existence of giant PBG [? ? ] and the tunable control of the band gap edges [? ? ] . 
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FIG. 1: Geometry of the problem. A Gaussian source is embedded in the center of an slab of high 
refractive index. The slab is bound with an arbitrary number of N periods at both sides. 

The control of fluorescence have been predicted on several PBG [? ]. In the present study, 
we investigate the change on the emitted power spectra due to the placement of a Lumines- 
cent dye within a flnite ID-PC. We have recently considered this problem theoretically to 
determine the degree of tuning of lasing due to the variation of the ambient condition in a 
liquid crystal based PC [? ]. The enhancement and suppression of the spontaneous emission 
have also been predicted and experimentally verifled in higher dimensional PBG structures 
[??]. 

In this work, we perform numerical modeling by using the Finite Difference Time-Domain 
method (FDTD) which is a methodology closer to the experimental conditions than the 
analytical formulations [? ? ]. We solve numerically the Maxwell equations using FDTD 
simulations performed by Meep, a freely available software package. [? ] This approach is 
flexible and also opens the possibility to account of defects or surface states [? ]. 

The geometry of the problem is shown in Fig. 1. A Gaussian source is embedded in a 
finite multilayer composed by alternative layers of high [uh) and low {11^) refractive indices 
illustrated in blue and red, respectively. The thickness of the slabs of high and low refractive 
indices are dn and d^. The period of the structure of the unit cell is d = dn + d^. The 
layers stacks above and below the source have an identical number of periods. The total 
thickness of the structure is L = 2Nd + dn- In the direction of the periodicity a detector is 
placed a distance s from the source. 
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FIG. 2: Relation dispersion and DOS in panels (a) and (b) for an infinite ID-PC of high and low 
refractive indices nn = 3.42 and til = 1.49, respectively. The width of the layers is dn = 0.4d and 
di = 0.6d, in each case. 

The dispersion relation of an infinite ID-PC is obtained using a well known formula [? ] 
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where kiuj) is the Bloch wave vector, and ku = uhoj/c and k^ = ulu/c are the wave 
vectors for the high and low dielectric media, respectively. It is expected that the radiation 
of the Luminescent dye embedded in the ID-PC can be related to the Density of States 
(DOS) which is defined as p{uj). p{uj) is obtained by differentiating the dispersion relation 
in the form [? ] 



dk{u) 
duj 



(2) 



The DOS is the number of allowed states and is obtained taking the relation between the 
number of modes Ak available for photons within the frequency Aw. 

A particular case of the relation dispersion and DOS is shown in Fig. 2 in panels (a) and 
(b), respectively. We take the values of the indices as rin = 3.42 and rii = 1.49. The width 
of the layers are dn = OAd and = 0.6d. We observe in panel (a) the existence of four 
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FIG. 3: SE from a Gaussian source embeeded in a finite ID-PC. The source is emitting at the 
center of the PBG with a central frequency ujc = 0.44 and a width of cr = 0.06. In panel (a) we 
present with solid line the emission for the case of the vacuum and for the eases of N=l,2,3 with 
dashed, dotted and tiny dotted, respectively. In panels (b) and (c) we present the cases of N=4,5,6 
and N=7,8,9. 



bands and three PBG. It has been reported by Joannopoulous et al. that the bands can 
be distinguished by where the power of the modes lies [? ]. The energy of the bands can 
be concentrate in the high index media or in the low index media. For this reason we refer 
to " high dielectric band" and " low dielectric bands" . A similar situation exist in electronic 
band structure of semiconductors where arc defined the "conduction" and "valence" bands [? 
] . We present in blue color the first and the third bands which are related to electromagnetic 
modes located in the the high index media. In the same manner, the second and the fourth 
bands are related to a electromagnetic mode located in the low index media. In panel (b) we 
present the DOS for each band that is obtained from the derivative of the relation dispersion 
using eq. (2). We observe an high value of the DOS at all the band edges. For clarity, we 
present in gray color all the PBG regions. 

The SE a Luminescent dye located in the center of the finite ID-PG is modeled as a 
Gaussian source for the electric field [Eg) in the form [? ] 
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FIG. 4: SE from a Gaussian source embeeded in a finite ID-PC. The source is emitting at the 
center of the PBG with a central frequency ujc = 0.48 and a width of cr = 0.06. In panel (a) we 
present with solid line the emission for the case of the vacuum and for the cases of N=l,2,3 with 
dashed, dotted and tiny dotted, respectively. In panels (b) and (c) we present the cases of N=4,5,6 
and N=7,8,9. 



Es(uj) = £;oe-°-5[<^-'-=H' (3) 

where cUc and a are the central frccuency and width of the pulse, respectively. 

The SE of a Gaussian source embedded in a Finite ID-PC are displayed in Fig. 3. We 
consider that the source have a central frequency of Uc — 0.44 and a — 0.06, respectively. 
In panel (a) we present as a reference the case of the emission in the vacuum with a solid 
line. We also present the cases of = 1, 2, 3. In panels (b) and (c) we present the emission 
for N = 4, 5, 6 and N = 7, 8, 9. We observe that as we increases the number of layers, the 
SE decreases. This is the expected results because the source is emitting in the PBG. 

In Fig. 4 we present the cases of the SE at the high energy limit of the band gap. We take 
a central frequency and pulse width of cUc — 0.48 and a — 0.06, respectively. We observe in 
panel (a) that the SE for the case of the vacuum and = 1 is almost similar, but it exist 
a dramatic enhancement of the radiation for N — 2,3. This is result of a redistribution 
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FIG. 5: SE from a Gaussian source embeeded in a finite ID-PC. The source is emitting at the 
center of the PBG with a central frequency ujc = 0.41 and a width of cr = 0.06. In panel (a) we 
present with solid line the emission for the case of the vacuum and for the cases of N=l,2,3 with 
dashed, dotted and tiny dotted, respectively. In panels (b) and (c) we present the cases of N=4,5,6 
and N=7,8,9. 

of the electromagnetic radiation induced by the existence of allowed Bloch modes. This 
preferential distribution of the electromagnetic field in the periodicity direction is similar in 
panels (b)and (c). 

The SE emission at the low energy limit of the band gap is presented in Fig. 5. Here 
we consider the values Uc = 0.41 and a = 0.06, respectively. In difference to the case of the 
high energy limit, here we do not observe an enhancement of the radiation measured by the 
detector, as it could be expected by the presence of the band gap edge. We observe in panel 
(a) a decreasing in the radiation. In panels (b) and (c) we observe a kind of stabilization of 
the radiated power. 

A global vision of the emitted radiation for the tree cases can be observed in Fig. 6. For 

the cases of the emission at the middle of the PBG {uc = 0.44) and at the lower band edge 
{uc = 0.41), we observe an important decreasing of the emission. In contrast, at the upper 
PBG {uJc — 0.48) we observe a high enhancement of the radiation. However, it is strange 
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FIG. 6: Comparison of the Emission as a function of the number of layers N for the cases cjc = 0.41, 
LOc = 0.44 and ujc = 0.48 displayed with red dotted, green tiny dotted and black solid, respectively. 

the low coupling of the electromagnetic field at the low energy limit of the band gap. Why 
the electromagnetic field is not coupled with this band? 

In order to explain this difference, we show in Fig. 7 the electromagnetic field produced 
at the low (wc = 0.41) and high {uc = 0.48) band edges in panels (a) and (b), respectively. 
We observe that in panel (a) exist a complicate mode, where a destructive interference exist 
that does not allow the excitation of the low index mode. In contrast, in panel (b) we 
observe that the high index mode is excited. The reason of this difference is that the source 
of electromagnetic field is placed at the high index material. 

In conclusion, we have presented a numerical modeling of the emission of a Luminescent 
dye embedded in a finite ID-PC. The treatment allows us to understand that is possible 
the enhancement of the radiation at the high energy limit of the band gap. We have also 
demonstrated that the DOS can not be directly related to an enhancement of the emitted 
power because it is necessary to take into account the conditions of coupling between the 
radiation of the source to the allowed Bloch modes at the band edge. 
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FIG. 7: Electromagnetic field at the band gap edges. In panel (a) and (b) we have the field 
distribution for the low and high limits of the band gaps, respectively. 
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